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Abstract 

The differential energy spectrum of cosmic rays from Cherenkov radiation measurements in EAS in the 
energy range of 10^^ — 10^*^ eV has been compared with an anomalous diffusion model for the particles 
in interstellar space having fractal properties (Lagutin et al., 2001). The close association between experi- 
mental data and calculated "all particle" spectra in form at Eq ~ (10^^ — 10^^) eV is found. In this case, 
the average mass composition of cosmic rays calculated by five components does not contradict the average 
mass composition from experimental data which was obtained by several of EAS characteristics in that 
energy region. 

I. INTRODUCTION ated directly with processes in interstellar space, 

namely, with the origin, acceleration and prop- 



agation of cosmic rays in the Galaxy and be- 
yond. The interpretation of these experimental 



The discovery of irregularities in the cos 

mic ray energy spectrum at the energy of ^ 

o imts 17 /T^i • .1 ^c^cr irTTiN J f^cts usiug thc diffcrcut models of cosmic ray 

3 X 10^^ eV (Khristiansen et al., 1956 yj) and ^ ^ 
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the detection of sharp decreases in the cosmic 
ray intensity at Eq > 5 x 10^^ eV (Greisen- 
Zatsepin-Kuzmin effect, 1966 (EHH) at the EAS 
arrays in Yakutsk, HiRes (USA), AUGER (Ar- 
gentina) are the most important achievement 

in the investigation of the superhigh and ultra- In this paper the comparison of the cos- 
high energy cosmic rays in recent years. Such mic ray energy spectrum by EAS Cherenkov 
a character of spectrum turn out to be associ- light measurements at the Yakutsk array E [U 

with the calculations according to an anomalous 

diffusion model of cosmic rays in interstellar 
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II. METHOD TO CONSTRUCT THE EAS 
SPECTRUM 

The showers at the Yakutsk array are selected 
with the central register by both scintillation and 
Cherenkov "masters" [fTOl fTTI . The all showers 
registered form the database of the Yakutsk EAS 
array. 

To construct the spectrum in energy, scat- 
tered by EAS particles in the atmosphere 
(Cherenkov radiation) the following selection 
criteria of showers are used: a) a shower core is 
to be placed within a perimeter of the array for 
the giant showers and near a center of the array 
for the showers with Eq < 10^^ eV. The show- 
ers whose cores are near the observation station 
i? < 60 m are excluded from a sampling: b) the 
probability to register a shower by Cherenkov 
photons is Wch > 0.9; c) a zenith angle is less 
than one-half of an aperture of Cherenkov de- 
tector, i.e. 9 < 55° in the case of the detector of 
the first type and 9 < 60° for the second type de- 
tector; d) a transmission coefficient of the atmo- 
sphere is > 0.60 for the wave length of 430 nm. 

Thus, more than 60000 showers with Eq > 
10^^ eV and 300000 showers with 10^^ < < 
10^^ eV were recorded in the database. To con- 
struct the spectrum, the showers were selected 
by the classification parameter Q{R = 150), 
i.e. by Cherenkov light flux density at a dis- 
tance 150 m from a core, which was propor- 
tional to the primary shower energy. The mea- 



surement accuracy for Q{R = 150) in the indi- 
vidual showers was 6 = AQ{R = 150) /Q{R = 
150) > 15%. 

The estimation of the shower energy Eq 
is determined by a quasicalorimetric method 
which does not depend on the EAS development 
model. A basis of the method is experimen- 
tal data about the Cherenkov light total flux, F, 
the total number of charged particles, N^^, and 
the total number of muons with i?thr > 1 GeV, 
A^M [[nilllini. The energy of individual show- 
ers is determined by the following formula: 

Eo = (9.1 ± 2.2) ■ 10^^ X Q{R = i50)0-99±o.o2 

(1) 

The intensity of cosmic ray flux in the 
given interval of EAS classification parameter 
is found as a ratio of the number of registered 
EAS events to S'eff ■ T ■ n. 

III. RESULTS AND DISCUSSION 

The differential energy spectrum of primary 
cosmic rays in the interval of 10^^ — 5 x 10^^ eV 
obtained from a totality of the all Cherenkov 
detector measurement data at the Yakutsk EAS 
array is shown on the Fig. [T] Our data con- 
firm an irregularity of the spectrum of "knee" 
type in the energy range of (2 — 5) x 10^^ eV 
discovered in [[U, and the irregularity of "an- 
kle" type at ~ 8 X 10^^ eV. It is estab- 
hshed that in the first case the spectrum in- 



dex is 7 = 2.7±0.1 below the break and 
7 = 3.03 ± OmatEo > 3 x 10^^ eV, and in 
the second case, the more sloping spectrum with 
7 = 2.6 ± 0.3 at > 8 X 10^*^ eV is observed. 
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Figure 1 : Energy spectrum of primary cosmic rays 
by measurement data of Cherenkov light at the 
Yakutsk complex EAS array. 

For the period of continuous observations of 
Cherenkov radiation more than 30 years (10% 
relative to one year time of EAS registration 
with the scintillation detectors), the showers 
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Figure 2: Differential cosmic ray intensity versus the 
energy. The points are Yakutsk array data, curves are 
the calculation from [9J. 



with Eo> 6x 10^^ eV did not detect. This fact 
confirms once more the GZK hypothesis flSl lU 
about the sharp break in the cosmic ray energy 
spectrum at .Eo > 5 x 10^^ eV. 

The galactic model 

The attempt to explain a form of obtained 
spectrum from the point of view of cosmic ray 
anomalous diffusion model and fractality of the 
Galaxy's magnetic field was made by Lagutin 
et al flU. The basis for the cosmic ray propa- 
gation in the Galaxy is the following assump- 
tions: a) after the generation in the sources, 
the particles move in fractal interstellar medium 
by means of two ways: the first way is "Levy 
flights", the second way is the motion along a 
spiral in the nonhomogeneous magnetic field, b) 
the particles exist during anomalous long time. 
The lifetime of particles is of a wide distribu- 
tion and its tail is described by a power law 
q{t) oc B ■ t-^-^,t ^ oo,/3 < 1 (so-called 
"Levy trapping time"). Calculations of the spec- 
trum were separately made for each of follow- 
ing groups of nuclei: p. He, CNO, N-Si, Fe. 
The resulting sum spectrum for the all particles 
is shown by a solid curve in Fig. [TJ From the 
calculations it follows that the suggested model 
reproduces the irregularity in the energy spec- 
trum of the "knee" type at Eq 3 x 10^^ eV 
and also the irregularity of the "ankle" type at 
i?o — 8 X 10^^ eV. This model does not ex- 



plain the behavior of a spectrum in the energy 
region of 10^^ — 10^^ eV and the break of the 
spectrum at Eq > 6 x 10^^ eV in more detail. 
The mass composition in the energy region of 
5 X 10^^ — 5 X 10^*^ eV is some heavier than at 
Eq ^ 10^^ eV, but this change is not very signif- 
icant that is expected from an experiment (see 
Fig. [3]). 

The galactic model with the sources of two 
types 

In contrast to [9], in the paper [fT5l a sce- 
nario is considered, in which supemovae are 
the main sources of cosmic rays and the ac- 
celeration up to -Emax — 105 ■ Z GeV takes 
place in the shock fronts. The particle spectrum 
formed in this case can be presented in the form 
of ~ E^'^6{E^a^ — E), where the Heavi- 
side function 9(x) reflects qualitatively the pres- 
ence of a sharp cut-off in the spectrum at > 
Emax [fT6l [TTll . The new calculations fulfilled 
by the above scenario of particle generation in 
the sources of two different types under the as- 
sumption of anomalous diffusion of particles 
in inhomogeneous medium show that at some 
parameters the anomalous diffusion model de- 
scribes satisfactorily the features of cosmic ray 
energy spectrum and mass composition up to 
Eq ~ 10^^ eV observed in an experiment. First 
of all, it refers to the fine structure of cosmic 
ray intensity change depending on the energy 
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Figure 3: Mass composition of cosmic rays at super- 
high and ultrahigh energies. The curve is a calcula- 
tion by Lagutin et al (2001) according to the anoma- 
lous diffusion model for cosmic ray propogation. 

(see Fig. [1]). By using these calculations, the 
sharp peaks in the mass composition depending 
on energy are also observed (see Fig. |3]). In this 
connection, it is of interest to compare calcu- 
lations in mass composition with experimental 
data obtained at the Yakutsk EAS big and small 
Cherenkov sub-arrays in recent years. 

Mass composition of primary cosmic rays 

Fig. [3] presents the results in mass composi- 
tion of primary cosmic rays of the Yakutsk ar- 
ray. The data were obtained in the framework 
of QGSJET-01 model and two-component mass 
composition (proton-iron nucleus). The several 
characteristics corresponding to the radial and 
longitudinal development of EAS are used in the 
analysis CH [Il|20l IZD 121 . 

The value (In A) in each case is determined 
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Figure 4: Comparison of the experimental spectrum 
with the calculated spectrum from ll23l for the meta- 
galactic protons (Eq > 5 x 10^^ eV) and galactic 
iron nuclei {Eo = 10^^ - 5 x 10^^ eV). 



by using the interpolation method [|24||. It is 
seen from Fig. [3] that the mass composition is 
varied up to heavy elements in the energy re- 
gion of (2 — 5) X 10^^ eV and becomes more 
light beginning with Eq 3 x 10^*^ eV. 

The lines are the calculations according to 
the anomalous diffusion model for the propaga- 
tion of cosmic rays in the Galaxy (Fig. [3} two 
sources) in the case of inhomogeneous galactic 
medium. In the first case, the monotone change 
in the mass composition up to -Eq > 3 x 10^^ eV 
is observed, after of which the mass composi- 
tion becomes more light. In the second case, the 
complicated structure in the dependence of mass 
composition on the energy is observed, peaks 
for the nuclei of different mass in the energy re- 
gion of 10^^ — 10^^ eV are noticeable. Accord- 
ing to a hypothesis [12511 and calculations from 
ifTSll . such an inhomogeneous structure can be 
formed by a near supernova. Our data (Fig. [3]) 



do not contradict to this hypothesis. 

Such a sharp change of the mass composition 
in the energy region of 5 x 10^^ — 5 x 10^^ eV 
is not explained in the framework of the galactic 
model and is likely associated with the existence 
of a transition boundary from galactic to meta- 
galactic cosmic rays. This conclusion is con- 
firmed by calculations from [|23l . where a sce- 
nario of galactic and metagalactic origin of cos- 
mic rays is considered. These calculations are 
shown in Fig. 4 together with our experimental 
data. It can be seen from Fig. |4]that cosmic rays 
in the energy region of 5 x 10^^ — 5 x 10^^ eV are 
most likely of galactic origin with the noticeable 
portion of heavy nuclei in the total flux. 

It should be noted the estimations of cos- 
mic ray mass composition in the region after the 
"knee", obtained at the compact arrays, agree 
well with each other. The same cannot be said 
of the energy region of ~ 10^® eV (see Fig. [s]) 
where HiRes array data point to more speedy 
enrichment of primary radiation by the light nu- 
clei and protons as compared with the Yakutsk 
array data. The Yakutsk FAS array data, on 
the contrary, show the gradual change from the 
heavy to light composition (protons and He nu- 
clei) at Eq ~ 10^^ eV. In both cases, data 
point to the existence of the tendency of "pro- 
tonization" of primary cosmic rays at Eq > 
3 X 10^^ eV. 



Conclusions 

Direct measurements of the cosmic ray en- 
ergy spectrum in the region of ultrahigh energies 
(in energy scattered by EAS particles in the at- 
mosphere) have confirmed the complicate form 
of spectrum. The spectrum becomes steeper 
at Eq > 3 X 10^^ eV and more sloping at 
.EO > 8 X 10^*^ eV. A character of energy de- 
pendence of {In A) by the Yakutsk EAS data 
point to the change of the mass composition 
of primary particles at singular points of cos- 
mic ray energy spectrum. The value (In A) rises 
with the energy after the knee up to its maxi- 
mum equal to 3.5 at (2 — 5) x 10^^ eV and then 
it begins to decrease. Such an energy depen- 
dence of (In A) does not contradict a hypothe- 
sis of cosmic rays propagation according to laws 
of the anomalous diffusion model in fractal in- 
terstellar medium (Lagutin et al., 2001). The 
value {\nA) at Eq > 10^^ eV decreases gradu- 
ally and at Eq ~ 10^^ eV the mass composition 
consists of He nuclei and protons. The cosmic 
ray intensity beyond Ethr. > 6 x 10^^ eV de- 



creases sharply and this effect is not described 
in the framework of the galactic model only. 
Such a character of spectrum does not contra- 
dict to the calculations by Berezinsky et al [|23ll 
for the metagalactic model, in which the "an- 
kle", observed in the experiments on ultrahigh 
energy cosmic ray registration, can be produced 
by the proton component only arriving from 
the Metagalaxy. Thereby, the details of exper- 
imental spectrum form, for example, "dip", i.e. 
the decrease of intensity at Eq x 10^^ eV, are 
caused by, most likely, the interaction of extra- 
galactic protons with a relic radiation photons 
(P + 7CMB ^ p + e+ + e"). As a direct argu- 
ment of this hypothesis, the anisotropy can be 
used which is related to the origin and sources 
of cosmic rays. Based on data of [|26l l27l l28l . 
at Eq > 8 X 10^^ eV the weak correlation in the 
arrival directions of EAS with the Galaxy plane 
and the close correlation with the Supergalaxy 
plane are observed and that the quasars can be 
the possible sources of cosmic rays. 
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